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^^ ■ We consider free and self-interacting quantum scalar fields satisfying modified dispersion rela- 

^^ ' tions in the framework of Einstein-Aether theory. Using adiabatic regularization, we study the 

^^ , renormalization of the equation for the mean value of the field in the self-interacting case, and 

^ . the renormalization of the semiclassical Einstein- Aether equations for free fields. In both cases we 

Q ' consider Bianchi type I background spacetimes. Contrary to what happens for free fields in flat 

•"T' ^ Robertson- Walker spacetimes, the self-interaction and/or the anisotropy produce non-purely geo- 

' . metric terms in the adiabatic expansion, i.e terms that involve both the metric g^^^ and the aether 

OO ' field u^. We argue that, in a general spacetime, the renormalization of the theory would involve new 

Cn I counterterms constructed with y^^ and it^, generating a fine-tuning problem for the Einstein- Aether 

theory. 
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I. INTRODUCTION 



In the last years it has been realized that the (still unknown) physics at very high energies may not be 



^^ ' inaccessible from an observational point of view. Indeed, trans-Planckian physics may have left an imprint 

If^ \ in the inhomogcneities of the cosmic microwave background radiation [l| , in the evolution of the scale factor 

lO ■ of the universe 0] , in the propagation of gamma ray bursts 0] , etc. 

^^ [ In the absence of a full theory, the theoretical approach to this problem is phenomenological. One possi- 

bility, that we will consider here, is to assume that the physics at high energies is such that its main effect is 
a modification of the dispersion relation of the quantum fields, thus violating Lorentz symmetry. Although 
t^^ , this is a simplistic approach, it could be useful to investigate whether the trans-Planckian effects could lead 

^^ ' to observable consequences or not in a given particular situation, by testing the robustness of the results 

under changes in the dispersion relations at very high energies. 
. , ■ The Modified Dispersion Relations (MDR) will obviously affect the structure of the quantum field theory, 

r> I in particular its renormalizability. Having in mind applications to cosmology, in previous papers [J, la, ly] , 

S , we have analyzed in detail the renormalization of free field theories with MDR in flat Robertson Walker 

spacetimes. We have shown that the theory can be renormalized using a generalization of the well known 
adiabatic regularization 0, Q that is used in theories with standard dispersion relations. As for the usual 
case, the adiabatic expansion of the energy momentum tensor contains divergent terms that can be written 
in terms of geometric tensors in n-dimcnsions, and therefore the theory can be renormalized by absorbing 
the infinities into the bare gravitational constants of the theory. It is remarkable that this can be done 
whatever the dispersion relation. This somewhat surprising result could be a peculiarity of flat Robertson 
Walker metrics [9[ and/or valid only for free flclds, and therefore it is of interest to investigate more general 
situations. 

In this paper we extend the adiabatic regularization to the case of self-interacting fields and anisotropic 
metrics (Bianchi type I). We will work within the context of the so called Einstein- Aether theory [lO|, a 
covariant theory of gravity in which the metric is coupled to a dynamical vector field. This field, that breaks 
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Lorentz invariance dynamically, is also coupled to the derivatives of the quantum matter fields, leading to 
MDR that contain higher powers of the momenta. The specific model is introduced in Section If. 

In Section III we consider a self-interacting scalar field on Bianchi type I metrics and discuss the renor- 
malization of the equation for the mean value of the field (f>Q. In order to do this, it will be necessary to 
compute the mean value of the fluctuations of the field {(fP) = {{(p — (f>o)^). We will calculate explicitly this 
quantity up to the second adiabatic order and show that, contrary to what happens for the usual dispersion 
relation, the second adiabatic order cannot be entirely written in terms of the metric and its derivatives, but 
also involve the aether field w^ and its derivatives. This property of {(f)'^) is valid even for free fields in flat 
Robertson- Walker spacetimes. 

In Section IV we analyze the renormalizability of the Semiclassical Einstein- Aether Equations (SEAE) for 
the case of free scalar fields with MDR in Bianchi type I universes. We compute (T^i/) up to the second 
adiabatic order. The zeroth adiabatic order is divergent whatever the dispersion relation. Being proportional 
to g/iu, the divergence can be absorbed into a redefinition of the comological constant. The second adiabatic 
order is shown to be divergent for dispersion relations that involve powers of the momenta smaller than 
or equal to four. This adiabatic order contains a term proportional to G^i^, that renormalizes Newton's 
constant. However, it also contains an additional non-purely geometric term, proportional to the variation 
of (V^u^)^. When this term is divergent, a new counterterm has to be introduced to renormalize the theory, 
even if originally not present in the classical Lagrangian. On the other hand, if it is finite, a counterterm 
would be necessary to make the theory consistent with observations. 

In Section V we argue that, for a general metric, the renormalization of the infinities produced by a 
quantum free field satisfying MDR will induce all possible counterterms involving up to two derivatives of 
the metric g^i, and the vector m^. As shown in Ref. [ll|j the coefficients of terms like (V^w'')^, R^uU^u" , 
etc, are strongly constrained observationally by post-Newtonian parameters, and therefore the counterterms 
induced by trans-Planckian physics should be fine tuned to satisfy these constraints. 

Throughout the paper we set c = 1 and adopt the sign convention denoted (+H — h) by Misner, Thorne, 
and Wheeler 111. 



II. THE MODEL 

We work in the frame of a generally covariant theory of gravity coupled to a dynamical vector field u^ 
that breaks local Lorentz symmetry. The most general action that is quadratic in derivatives is given by 

Sg = -^ I rxV^iR - 2A + £„), (1) 

where g = det{gfj_,y), R is the Ricci scalar, A and G are the bare cosmological and Newton's constants, and 
£u describe the dynamics of the additional degree of freedom u^, 

Cu = -ACg^^u^u, + 1) - fei^^.F^" - 62(V^u^)2 - b^R^.u^u^' - h^uPu^V pU^V ^u^', (2) 

where F^^i, = V^Ui/ — V,yU^. The Lagrange multiplier A is introduced to impose the condition u^u^^ = —1 
and the coefficients hi {i = 1, 2, 3, 4) are arbitrary. The term V^w^Vw^ coincides with (V^m^)'^ — R^^u^u'^ 
up to a total derivative, and hence has been omitted. 

We consider a quantum scalar field (j) with a generalized dispersion relation propagating in a curved 
space-time with a classical background metric given by 

ds^ = gf.^dx'^dx'' = -{u^dxf'f+ L^^ dx^^dx", (3) 

where fi,p = 0, l...n — 1 (with n the space-time dimension) and -Lf^i,= g^,^ + u^u^. The action for the scalar 
field can be written as: 

5*0 = / d^Xyf^iC^ + Ccor + C,nt), (4) 



where C^ is the standard Lagrangian of a free, massive, minimally coupled scalar field 
Ccor is the corrective lagrangian that gives rise to a generalized dispersion relation 



(5) 



(6) 



where < p < s, hsp are arbitrary coefficients, and T) (j) — J-u ^a -L^ V'^ (V^ is the covariant derivative 
corresponding to the metric g^i, and -Lu= 9^'^ -L/ji^)- The interaction Lagrangian Cmt contains the following 
terms: 

where ^, S.i {i = 1, 2, 3, 4, 5) and A are bare parameters. Note that, in addition to the self-interaction and the 
standard coupling to the Ricci scalar, we have also included couplings between t/)^ and non-purely geometric 
terms that involve the aether field u^. Note also that, if we assume that the MDR depart from the usual 

one at a given scale Mc^ the coefficients bgp scale as bsp ^^ M(j ^ ^ . 

In the rest of the paper we will consider a four-dimensional Bianchi type I space-time with line element 



ds^ = -dt^ + Y^ C,{t)dxl = ~C{ri)df]^ + ^ C^{t)d2 



(8) 



i=l 



where C = (CiCaCa)!/^^ di] = dt/C^/'^, and u^ = C'^/^{'n)5l. Therefore, in this frame F^^ = and 
w^V^u^ = 0. In what follows we use primes for denoting derivatives with respect to the conformal time rj. 
No sum convention in spatial (latin) indices is assumed. The generalized dispersion relation takes the form 



C(.?) 



+ x + 2Y^{-~\Y+Pbspx'^'+P^ 



(9) 



where x = Eti fef/Q = ELi ^^ ^ ELi ^Af , with Eti Af = 1. 



III. SELF-INTERACTING SCALAR FIELD IN BIANCHI TYPE I SPACE-TIMES 



In this section we are concerned with the rcnormalization of the equation of motion for the expectation 
value of a self-interacting scalar field (A 7^ 0) propagating in a four-dimensional Bianchi type I space-time. 
We assume that the state of the system is such that the expectation value of the field is 0o . Then, defining 
a new quantum field (j) as <() = 4)q + 4>, the equation of motion for (f>Q in the one-loop approximation is given 

by 



□0 



s,p<s 



The Fourier modes of the scaled field x — C^^^(i> satisfy 
XfcH 



1+ U-lj CR + Q + 6C^(V^M^)' + ^sCV^u^V'u^ + i^CR^^.u^u" + UCXcg 



00 - 4A03 ^ 0. 
(10) 

Xk = 0, (11) 



with the usual normahzation condition 



XkX'k* - x'kXl = i 



The exphcit expressions for the different terms in Eqs. pO|) and pT|) are, in Bianchi type I metrics, 
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where d; = C[/Ci and D = X]i=i '^i/3 = C jC. Note that for the metric we are considering 



(12) 

(13a) 
(13b) 

(13c) 
(13d) 

(13e) 
(14) 



and therefore without loss of generality we can set ^5 = 0. 

For dispersion relations such that the mean value (0^) in Eq. ([rU]) is divergent, the infinities must be 
absorbed into the bare constants of the theory. To implement the renormalization, we start by expressing 
the field modes Xk in the well known form 



1 



Xk = 



V2Wj 



exp —I 



Wkif])dfj) , 



which allows us to write 



1 



'^ ' (27r)3C 
Substitution of Eq. US]) into Eq. ITT]) yields 
1 



d^k\xk\^ = 



(27r)3C 
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(15) 



(16) 



(17) 



For adiabatic regularization we need the approximate solution of this non-linear differential equation that is 
obtained by assuming that W^ is a slowly varying function of 77. In this adiabatic or WKB approximation 
the adiabatic order of a term is given by the number of time derivatives of the metric plus the power of (J)q 
[13| . The WKB approximation can be obtained by solving the Eq.([T7]) iteratively 
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(18) 



where the superscript denote the adiabatic order. To lowest order we have ^"^VF/t = ujk- The second adiabatic 
order can be computed replacing Wk by tok on the right-hand side of Eq. (fTT)) . Thus, we straightforwardly 
obtain 



^^^Wl = Ci?(e-i)+Q + ^-^+6C^(VX)'+e3CV^w.V^t^ + 12AC0^ 
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where we have defined the function 
with ljI = ujI/C. We have also used that 

,2^/ 3 
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dlnCj 
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(21a) 
(21b) 
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^j=l 



where a dot indicates a derivative with respect to In a;. 

We proceed as for the standard dispersion relation, defining the rcnormalized expectation value as 



/ad2, 



(22) 



with U^)ad2 = ((/'2)(0) 



j)'^)^'^\ where again the superscripts indicate the adiabatic order. 

We now compute the zeroth adiabatic order of {(j)^) and regularize it by using the fact that the integral of a 

total derivative vanishes in dimensional regularization \XM- ^'^^ this, and in order to avoid the complications 

of computing all quantities in n-dimensions, we first perform the angular integrations and then generalize 

the four-dimensional integrals to n-dimcnsions by replacing d^k = C^^^d^y = C^/^y^dydVt {yi = ki/y/Ui) by 

Therefore, the zeroth adiabatic order is given by 
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where /i is given Table ID Note that the integral /i is divergent unless uxf. behaves as x^ with s > 3, for large 
values of x. This divergence can be absorbed in the bare mass of the quantum field (see below). 
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TABLE I: Explicit expressions for li. To obtain these integrals we have made the change of variables x = y and we 
have defined ujk ~ LOhjyfC- 



The second adiabatic order can be written as 

/^2\(2) ^_ 
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The angular integrations can be performed with the use of the identities hsted in the Appendix A. After 
some calculations we obtain: 
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where I3 is given in Table [H and we have defined the integrals 
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with m, n, I, s integer numbers. As it is shown in the Appendix of Ref. [8 
terms of the ones in Table [J by performing integrations by parts. For n 

^1000 =0, 
2 



Then, substituting this results into Eq. 



■^2000 
•^0100 
3 we arrive at 
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this integrals can be expressed in 
+ 4, we have Q: 
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This is the main result of this section. The relevant point is that, in addition of the usual terms proportional 
to R and (J>q, the second adiabatic order contains terms with two derivatives of the aether field, which are 
present even if ^2 = Cs — 0. For the standard dispersion relation I3 diverges and I4 vanishes (sec Table |T|. 
Therefore, when ^2 = ^3 = one reobtains the usual result. However, for any other dispersion relation of 
the type given in Eq. ©, I3 and /4 are finite. An interesting point is that, if we consider a generahzed 
dispersion relation, evaluate the integral explicitly in four dimensions and then take the limit in which the 
dispersion relation tends to the usual one, a nonvanishing finite result can be obtained. For example, a 
dispersion relation of the form uif = C{x + 2biix^) yields 



2bn 



+ 00 



dx{l + 2biixy2 



(29) 



Therefore, there is a finite remnant of the trans-Planckian physics in the second adiabatic order, even in the 
limit in which the scale of new physics is very high Mc -^ 00 {bn -^ 0). 

Coming back to the mean value equation ()10p . we write the bare parameters in terms of the renormalized 
ones plus the corresponding to countcrtcrms: 



acPo - [ml + 6m^ + {^r + S^R + {^2R + <56)(V^u^)2 + (^3^ + 6^3)^ ^u^S/^u^ 
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00 - 4(Ai^ + SXm = 0. 



(30) 



Introducing Eqs. (^5]) and (^5]) into Eq. ([5(I)) . we see that the regularized second adiabatic order ((j)'^)ad2 can 
be absorbed into the bare constants by defining counterterms such that 

q\2 
^A=^/3, (31b) 

S^ =#^ fCi. - ^"l /3, (31c) 



(2^)2 V 6 

^6 =^^3 - ^/4. (31e) 

Note that, even when the parameters ^2i?, and ^3^ are set to zero, the corresponding counterterms arise due 
to the self-interaction of the scalar field. Note also that by considering the same theory but in a background 
flat FRW space-time, it is not possible to distinguish between the redefinitions of ^2 and ^3 proportional to 
/4, since in such background we have that 

W^u,W''u^ = ^{W^un^. (32) 

Prom the results of this section we conclude that, as long as one considers the renormalization of the 
mean value equation in Bianchi type I spacetimes, and for the class of MDR considered here, it is enough to 
subtract the zeroth adiabatic order of ((/)^), since the second adiabatic order produce a finite renormalization 
of the bare constants of the theory. It would be interesting to check whether this is a general property, i.e. 
valid for an arbitrary background, or not. In order to address this issue, it would be necessary to know the 
singularity structure of the two-point function of a quantum field satisfying MDR for arbitrary values of g^^ 
and w^. This singularity structure could be revealed by a generalized momentum-space representation of 
the Green's functions [la, [lg|. In any case, the calculation of the second adiabatic order presented in this 
section shows that the interaction terms proportional to ^2 and ^3 that appear in Eq. ([7]) are generated by 
quantum effects, even if not present at the classical level. It is likely that the other interaction terms will 
also be generated in a more general background. 

IV. ON THE RENORMALIZATION OF THE STRESS TENSOR IN BIANCHI TYPE I 

SPACE-TIMES 

In this section we focus on the renormalization of the SEAE. We restrict the analysis to the case of a free 
scalar field (A ~ 0, {(/)) — 0) and, for the sake of simplicity, we set the parameters (,i = (1=1,2,3,4,5). 
The SEAE take the form 

G^. + Ag^. = 8^G[t;,^ -f {t;:^ + T^t) + t;;^ (33) 

where A and G are the bare cosmological and Newton's constants, G^,y is the Einstein tensor, T^u = 
- -2= ^S^ and T", = T"5 + T,),'= , T"' is the stress tensor of the background vector field while T^,", is 

the additional contribution due to the modification of the Lagrange multiplier A arising from the coupling 
between the scalar field (f) and u^. T/;^"" is a stress tensor coming from classical sources not coupled to the 
aether field. As we will compute the mean value of the stress tensor up to the second adiabatic order, we 
omit classical terms quadratic in the curvature (we will comment on this issue in the next section). 



The nontrivial components of the Einstem tensor are, in Bianchi type I spacetimes: 
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The stress tensor corresponding to the background vector field u^ can be written as 
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whose nonzero components are 
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The expectation value of the quantum energy momentum tensor T^,j = T^^^^ + T^^ is given by 
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Using the expression given in Eq. (jlSp for the modes Xk^ it can be written as 
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Therefore, the zeroth adiabatic order can be expressed in the form 
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where we have used that / dflXf = 47r/3. Then, after an integration by parts in Eq. (j41b[) we obtain, as 
n ^ 4, 

where /o is a divergent integral as 77, ^ 4 for any of the dispersion relations given in Eq. ^ (see Table |I| . 
Hence, this regularized adiabatic order can be absorbed into a redefinition of the bare cosmological constant 
A. 

The second adiabatic order of (T^i/) can be written as 



(,) _ [dmyy^-^U(i \[{u^ MY (ojIy f_D e 



where ^^''W^ is given by the expression in Eq. (fT9|) with A = ^2 = ^3 = 0. The explicit expressions for 
{ioly/ivl and (a'^)"/wfe are given in Eq. ((2T|l . 

After performing the angular integrations with the use of the identities given in the Appendix A and some 
algebraic manipulations, we obtain: 

(^w)^'^ = TTTI ["1^' + "2Q], (45a) 

The coefficients a^ and /3i are given in Appendix B, where it is also shown that using integration by parts 
they can be expressed in terms of two of the integrals in Table U] Thus, wc find 
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^■i^4)-S 



G,,. + ^t;, } . (46) 



Note that both /i and I2 diverge when w^, behaves as a;*,s < 3 for large values of x. Therefore, in this 
case the divergences should be absorbed into the bare constants G and 62- However, when s > 3, the second 
adiabatic order produce finite renormalizations of both constants. 

As in the evaluation of (0^) presented in the previous section, depending on the dispersion relation one 
could have a remnant of the trans-Planckian physics in the second adiabatic order of (T^u)- Indeed, while 
I2 vanishes for the standard dispersion relation, a non vanishing (and even divergent) result can be obtained 
for MDR in the limit Mc -^ 00. For example, for a dispersion relation of the form w^ = C{x + 2b22X^) we 
find that 

/•+°° x^ of 
I2 = 24622 / dx = ^r [1/6] r [4/3] , (47) 

7o (1 + 2522x3)2 ^q_^ 

which diverges as Mc — *■ 00 (622 — *■ 0). 

Eq. (|^S)) is the main result of this section. We sec that, for a generalized dispersion relation of the 
type given in Eq. ([9|), not only a redefinition of the Newton's constant is necessary in order to cancel the 
divergences of the second adiabatic order, but also a redefinition of the coefficient 62 which corresponds to 
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the term (V^it'^)^ in the bare Lagrangian of the vector field. The second adiabatic order contains terms 
that are non-purely geometric, in the sense that they cannot be written only in terms of the metric, but also 
involve the aether field. 

It is noteworthy that for a background fiat FRW space-time G^y ~ 3/27"" , thereby, in Refs.0, Q it was 
not possible to realize that a redefinition of the Newton's constant is not enough for cancelling the second 
adiabatic order. In fact, for this particular space-time G^n is the unique covariantly conserved tensor of 
adiabatic order two that can be derived from an action formed by combining the vector field u^ , the metric 
g^jy, and their derivatives. 

V. DISCUSSION 

In this paper wc have worked within the context of a generally covariant theory of gravitation coupled 
to a dynamical time-like Lorcntz-violating vector field. Wc considered a quantum scalar field satisfying 
MDR, and analyzed the renormalization of the infinities that arise in the semiclassical theory. In particular, 
considering Bianchi type I spacetimes, we have analyzed the dynamical equation for the expectation value 
of a self- interacting scalar field (Section III), and the SEAE for the metric in the case of a free scalar field 
(Section IV). With the use of adiabatic subtraction and dimensional regularization, we have shown that, in 
addition to the usual terms required to absorb the infinities of the second adiabatic orders, it is necessary 
to consider more general counterterms that involve the aether field. This property was not apparent in our 
previous works [J, Q, due to the high symmetry of the fiat Robertson Walker metrics. 

These results suggest that, in a more general background metric, any covariant term which can be formed 
by combining the vector field m^, the metric 3^1/ and up to two of their derivatives, will appear in the 
regularized second adiabatic order of the expectation value of the quantum stress tensor, provided that the 
theory contains a scalar field with a generalized dispersion relation of the type given in Eq. ^. Hence, in 
order to absorb the divergences contained in the second adiabatic order, a bare action as general as the one 
given in Eq. ([T]) should be considered. Depending on the particular dispersion relation of the quantum field, 
the second adiabatic order may be finite. If this is the case, quantum effects generate finite renormalizations 
of the constants appearing in the classical Lagrangian. As wc have also pointed out in Section IV, this finite 
renormalizations could be extremely large. 

In the weak-field limit, the terms proportional to the constants bi in Eq. ^ could have observable 
consequences. Indeed, the most general action given in Eq. ([1]) has four free parameters more than general 
relativity. This theory has been studied in several contexts, such as of the static weak-field limit [ll|, the 
radiation and propagation of the aether-gravitational waves [l3|, cosmology [l^, etc., in which stringent 
constraints on the parameters have been imposed to make the theory consistent with observation. For 
example, in Ref. [11| it is shown that for all the PPN parameters to agree with observation, the four 
additional parameters of the model must satisfy two constraint equations with sufficient accuracy (i.e., the 
additional four-parameter space of the model has to be practically reduced to a two-dimensional subspace). In 
the absence of a known mechanism to explain why the parameters satisfy precisely such constraint equations, 
it seems that quantum effects generate a fine-tuning problem in the Einstein- Aether theory. This is analogous 
to the fine-tuning problem present in the Myers- Pospelov modification of QED [I3 . 

In this paper we restricted ourselves to the evaluation of the adiabatic expansion up to the second adiabatic 
order, in a particular class of background metrics. By power counting, wc expect the fourth adiabatic order 
(T'pi/)^*-' to be finite in these metrics. However^ there could be some subtleties related to the would be 
Gauss-Bonnet invariant in four dimensions [a, Q. In the light of the results obtained in this paper, this 
issue should be reexamined. One should compute (T^^)^"*^ for a general background metric and aether 
field (this could be done by generalizing the momentum-space representation of the Green's functions). 
On dimensional grounds we expect {T^uY^ to contain terms proportional to the variation of i?(V^u'')^, 
{Rfj.^u'^u'')^, Rf^t^u'^u" R, Rf^^pcru'^uPR'"', etc., in addition to the usual ones: i?2, R^^Rt"" and i^^^p^i?''"'"". 
Depending on the MDR, the fourth adiabatic order could be finite or divergent when expressed in terms of 
such variations. This fact would define whether it is necessary or not to subtract the fourth adiabatic order 
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in a general background, for a given dispersion relation. Work in this direction is in progress. 

Appendix A: Identities for Bianchi type I space-times 

In this Appendix wc briefly summarize some useful formulas required for the adiabatic regularization of 
(i^^)*^^' and (Tpi/)''^-' in Bianchi type I space-times. 

As we have already mentioned in the text, in order to regularize the theory we perform the four- dimensional 
angular integrations and then generalize the integrals to n-dimensions. We rescale the integration variables 
ki —^ yt = ki/Ci and transform the volume element (Py from rectangular coordinates to spherical coordinates 
y^dydfl, where dfl is the solid angle element. In terms of yf = y^Xf, the relevant integrals are of the form 

I{i,j,k) = JdnXfxl^Xf, (48) 

which can be evaluated by using the fact that they are invariant under permutations of {i,j, k}. We provide 
here a list of the integrals we have used in this paper (see [20| for more details): 

An 
mO,fc)=^^^, (49a) 

7(1,1,0)=^, (49b) 

7(1,2,0)=^, (49c) 

These results, together with the formula J2i=i '^1 = 3(8Q + D^), allow us to derive the following identities: 



J2^ f dndjdkX^.Xl =47r (d^ + ^Q 

3 = 1 k=l-' ^ 

±J:ldnd,d,Xi)^^xl^^' 

A — 1 7 1 ^ 



3 3 _ ^^^ 

l2 I ,1 J. n I ° J2 



j=l k=l 



bD^ + U,D + -df + 16Q 



(50a) 
(50b) 



V / dn{d'^ + 2d,D - d])XlX] =-^ [2{d[ + 2d,D - df) + 3(D' + D^ _ gg)] ^ (50c) 

that are useful for the evaluation of (c^^)'^' and {T^y)'^'^\ 

Appendix B: Regularization of {T^^Y'^' in Bianchi type I space-times 

In this Appendix we provide some details for computing the second adiabatic order of the expectation 
value of the quantum energy momentum tensor. The explicit expressions for the coefficients appearing in 
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Eq. (|i5l) are: 

ai =^[-4/io + ho + 4(/oo - ee/oo + 6^/10)], (51a) 

a2-^[/oo + 2/io+/2o-30e/oo], (51b) 

/3i =^[52/00 - 120/01 - 76/io + 20/11 + 77/20 - 5/30 + 280e(-/oo + 2/oi + ho - ho)], (51c) 

/32 =7j^[2/oo + 39/10 - 8/20 - 20^(22/00 + ho)l (51d) 

640 

/33 =77^ [-8/00 + 12/01 - 19/10 + 12/11 - 14/20 - 3/30 + 210e/io], (51e) 

lOoU 

Pi =Yi^[- 24/01 + 3/10 + 4/11 + 21/20 - /30 - 14e(-8/oi + ho + 4/20) + /oo(42e - 19)], (51f) 

/35 =^[2% + 4/01 + 3/10 + 4/11 - /30], (51g) 

/Se =Y^[-^oo - 2/10 - ho + 3OC/00], (51h) 
where the integrals Imn are given by 

dx^rr, (52) 

with m,n = 0, 1,2,3. 

Let us now sketch the procedure to find relations between these coefScients in the context of dimensional 
regularization, which is completely analogous to the one described in the Appendix of Rcf. 6] for relating the 
integrals Jmnis of Eq. (|26p . By definition, ho = h , and 

ho = / dx— — -Tn-r- - 1 = -2 / dxx^^^ h 53) 

Jo ^k \^t dx J Jq dx 

= (--iKi- ^1(^^2/1, 

where we have performed an integration by parts and discarded the surface term. Similarly, one can prove 
that 

(54a) 

(54b) 

(54c) 

(54d) 

where the integrals li on the right hand side arc given in Table IH and / (which docs not appear in the final 
results) is given by 



r 
20 


^., 


30 


4-1^^ 


r 
01 


= -27, + i/,. 


'11 


2 2 
-15^^+15^ 



k 



X 2 d^Cb'^ 



d^^^-nr- (54e) 
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Replacing these results into Eq. ([5T|) . we obtain: 



"^4^ (^4)+ 1' ^''^^ 

/?5 =0. (55g) 

Finally, after substituting these coefficients into Eq. (|l5|) we arrive at Eq. (|46p . 
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